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ABSTRACT 


Context. The north-west photo-dissociation region (PDR) in the reflection nebula NGC 7023 displays a complex structure. Filament¬ 
like condensations at the edge of the cloud can be traced via the emission of the main cooling lines, offering a great opportunity to 
study the link between the morphology and energetics of these regions. 

Aims. We study the spatial variation of the far-infrared fine-structure lines of [CII] (158 pm), [OI] (63 and 145 pm). These lines trace 
the local gas conditions across the PDR. We also compare their emission to molecular tracers including rotational and rovibrational 
H 2 , and high-rotational lines of CO. 

Methods. We use observations from the Herschel/PACS instrument to map the spatial distribution of these fine-structure lines. The 
observed region covers a square area of about 110"xll0" with an angular resolution that varies from 4" to 11". We compare this 
emission to ground-based and Spitzer observations of H 2 , Herschel/SPIRE observations of CO lines, and Spitzer/IRAC 3.6 pm images 
that trace the emission of polycyclic aromatic hydrocarbons. We use a PDR code to model the [OI] 145/rm line, and infer the physical 
conditions in the region. 

Results. The [CII] (158 pm) and [OI] (63 and 145 pm) lines arise from the warm cloud surface where the PDR is located and 
the gas is warm, cooling the region. We find that although the relative contribution to the cooling budget over the observed region is 
dominated by [O I]63 ^m(>30%), H 2 contributes significantly in the PDR (~35%), as does [CII] 158 pm outside the PDR (30%). Other 
species contribute little to the cooling ([OI] 145 pm 9%, and CO 4%). Enhanced emission of these far-infrared atomic lines trace the 
presence of condensations, where high excitation CO rotational lines and dust emission in the submillimitre are also detected. The 
[OI] maps resolve these condensations into two structures, and show that the peak of [OI] is slightly displaced from the molecular 
H 2 emission. The size of these structures is around 8" (0.015 pc), and in surface cover about 9% of the PDR emission. We have 
tested whether the density profile and peak densities derived in previous studies to model the dust and molecular emission, can 
predict the [OI] 145 pm emission. We find that the model with a peak density of 10 6 cm 3 , and 2xl0 4 ~ 5 cm -3 in the oxygen emitting 
region, predicts an [OI] 145 pm line which is just 30% less than the observed emission. Finally, we have not detected emission from 
[NII] 122 pm, suggesting that the cavity is mostly filled with non-ionised gas. 

Key words. Infrared: general - photon-dominated region (PDR) - ISM: lines and bands - ISM: individual objects - NGC7023 


1. Introduction 

NGC 7023 is a bright and well-known reflection nebula. 
Although it is a relatively nearby object, its reported distance 
varies significantly in the lite rature. Using Hipp arcos parallaxes 
Ivan den Ancker etaT. ( 119971) and later on Ivan Leeuwenl ( 20071 ) 
measured a distance of 430 4 ( l )( f j° pc and 520+180 pc respectively, 
and iBenistv et al.l (1201 3l) estimated a distance of 320+51 pc 
based on revised orbital paramete rs and astrometry. In this pa- 
per we adopt the distance given by Ivan den Ancker et al.l (119971) 
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(430 pc), which is almo st the average of th e values given by 
Ivan Leeuwenl (120071) and IBenistv et alJ (120131) . 

It is excited by the Herbig Be binary star system HD 200775 
dAlecian et all 120081), which i s imm ersed in the cavity of the 
molecular cloud dFuente et all 119921) , The overall morphology 
of NGC 7023 is driven by the action of the ionising stars, where 
the stellar outflow, now inactive, creates a biconical cavity in 
the_shape of a_butterfly of about of 1.5 pcx 0.8 pc in size 
dKohler et all 12014 s ). The stars illuminate the walls of the cav¬ 
ity in an almost edge-on orientation. At the edge of these walls, 
three PDRs are present dFuente & Martm-Pintadol 119971) . and 
are located at 40" northwest (herein NW), 70" south-west (SW), 
and 170" east of the stars. The NW PDR is the brightest, and 
is oriented almost edge-on to the observer. Polarization ob¬ 
servations in the near-infrared dSellgren et all I1992), as well 
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as vibrational Hi (iLemaire et aU 1995) of the NW PDR re¬ 
veal a filamentary morphology !Witt et al.l (20061) find that this 
morphology is the result of several condensations that are su¬ 
perimposed along the line of sight, whose surfaces are di¬ 
rectly illuminated by the star. The incident radiation field, in 
units of the interstellar far-u ltraviolet radiation fi e ld estimated 
by H abin JI is G =2.6xl0 3 (IChokshi et all fl988t iPilleri et all 
i2012t) . However, higher values also having been reported in the 
literature dFuente et al.L 119991) . Deriving the incident radiation 
field relies on several assumptions on the geometry and spec¬ 
tral type of the star. For the rest of the paper, we will assu me 
the value of G =2.6xl0 3 at the PDR front as in Choksh i et alJ 
dl988h an d lPillerietal] (120121) . 


NGC 7023 has been the subject of many studies in the lit¬ 
erature, including several spectroscopic studies in the mid- and 
far-infrared (FIR). The mid- and far-IR spectral regions con¬ 
tain many important cooling lines and hold the fingerprints 
of the dust emission. iFuente et al.l d2000t) presented ISO/SWS 
and LWS observations of NGC 7023 and differentiated between 
three regions: 1) the star and the cavity formed by the star, which 
is filled with a low density gas, 2) the edges of the cavity that 
defines the PDRs, and 3) the molecular gas. Based on their anal¬ 
ysis of the atomic and molecular lines ([O I]63 gm, [OI] 145 gm, 
[CII]158/rm, H 2 ), they find that the NW and SW PDRs have 
similar excitation conditions, with filaments having a high den¬ 
sity of n ~10 6 cm -3 , and a lower density of n ~10 4 cm~ 3 in 
the region between f ilaments. IWerneTet al. (2004) and, later on, 
Sellg ren et al.l (120071) studied the variation of the dust features 
in the mid-IR spectrum of NGC 7023 using the Spitzer/IRS and 
MIPS instruments. They observed variations in the NW and SW 
PDRs of the intensity, widths, and central wavelength of spectral 
features arising from polycyclic aromatic hydrocarbons (PAHs), 
as a function of the distance to the ionising source (HD 200775). 
IWitt et al! (2006 ) found a good correlation between the opti¬ 
cal emission from the dust, as traced by the R band (658 nm), 
with that of ro-vibrational H 2 . Using Spitzer data, [Fleming et af] 
( 2010 ) studied the H 2 and PAHs emission in the PDRs of three 
reflection nebulae, including NGC 7023. They showed evidence 
for PAH de-hydrogenation, which they concluded is suggestive 
of Hi formation on PAHs. 

Recently, NGC 7023 has been the target of several pr ogram s 
using the Herschel Space Observatory (iPilbratt et al.L [201 ()t) , 
vastly improving the spatial resolution of previous far-infrared 
studies. lAbergel et al.l (20101) combined Spitzer and Herschel 
maps to study the spatial variations of the dust properties in 
the east PDR. Using a radiative transfer code they were able 
to reproduce the different spatial variations of PA Hs^ very small 
grains (VSGs), and large dust grain emission. iJoblin et al.l (20101) 
used the high spectral resolution provided by the HIFI instru¬ 
ment to study the [CII] 158/mi line at several positions in the 
nebula, including a central position in the NW PDR. They found 
that the emission of [CII]158/mi line and PAHs arises from 
the transitio n region between the atomic and molecular gas. 
lOkada et al.l ( 2013 ) studied the photoelectric efficiency from 
PAHs in a sample of six PDRs, including the East and NW 
PDRs of NGC7023. Using Spitzer and Herschel/PACS obser¬ 
vation, they evaluate the photo-electric effect from PAHs and 
the cooling [OI]63 and 145 //m, and [CII] 158 /mi lines in three 
regions (cavity, interface, and molecular region) and find that re¬ 
gions with high fraction of ionised PAHs have a lower heating 
efficiency. 


1 The Habing field correspond s to 1 .6x 10 3 erg s 1 cm 2 when inte¬ 

grated between 91.2 and 240 nm (ITabi n ul l l 968 ). 



Fig. 1. Overlay of the PACS field of view on a ground based H 2 
1-0 S(l) map at 2.12 //m (Lemaire et all ll996l) . A 3x3 raster 
map was performed (see Sect. 2). The map covers an area of 
approximately 110"xll0" and is centred in the NW PDR. The 
emission from the ionising stars (HD 200775), on the bottom left 
of the figure, is saturated. RA and Dec coordinates are in degrees 
(J2000) 


In this paper we present observations from the 
Photoconductor Array Camera and Spectrometer (PACS, 
iPoglitsch et all 1201 Oh of the [OI]63/mn, [OI] 145/mi, and 
[Cn]158/rm lines in the NW PDR of NGC 7023. These obser¬ 
vations allow us to resolve in detail the global shape of the PDR, 
although substructures su ch as the filaments obs e rved at higher 
spatia l resolution (e.g., ILemaire et al., 1996; IFuente et al. L 
IT9961) . are unresolved. This study complements a recent SPIRE 
spectroscopic study of the molecular gas and dust emission 
by our group (Kohler et all 12014 ). This paper is organised as 
follows: the observations and data reduction are described in 
Sect. 2. In Sect. 3 the spatial morphology of the lines, their 
comparison to the emission of CO and Hi, and the correlation 
between the different line ratios as a function of the distance to 
the ionising source, are discussed. In Sect. 4 the cooling budget, 
and the morphology of the main cooling lines are analysed. In 
Sect. 5 the results of the modelling are described. Finally, the 
conclusions are summarised in the final section. 


2. Observations and Data Reduction 

The observations were carried out during the science demonstra¬ 
tion phase (SDP) of Herschel on 7 lanuary 2010, and are part 
of the ”Evolution of the Interstellar Medium” guaranteed time 
key project (observation ID=1342191152. lAhergel et al.U2010l) . 
The observations were taken using the PACS instrument in the - 
now decommissioned - wavelength switching mode. Four fine- 
structure lines were observed: [CII] at 158/mi, [OI] at 63 and 
145//m, and [Nil] at 122/mi. The observations were centred on 
the NW PDR, with RA and DEC coordinates of 315.375°and 
60.178°respectively (J2000). 

The observing strategy and redu ction methods follow that o f 
our earlier paper on the Orion Bar ( Bernard-Salas et all 1201 2). 
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We thus summarise the main steps here, and refer to the above 
paper for details. 

To trace the NW PDR a 3x3 raster map was performed. The 
configuration at the time of the observation is shown in Fig. 1, 
where the raster map at the epoch of observation is overlaid on 
top of an H 2 map of the region made with the Canada-France- 
Hawaii Telescope (CFHT). An additional map centred in the 
East PDR (not shown in the figure) will be presented in a fu¬ 
ture paper. A minimum exposure configuration of one cycle and 
repetition per line was performed. Raster and point line steps of 
23.5" were chosen and result in Nyquist sampling for the lines 
in the red channel ([CII] 158/mi, [OI] 145/mi, [Nil] 122/tm). 
The purpose of the wavelength switching mode is to cancel 
out the background by determining a differential line profile. 
This means that observations at an off-position are not needed 
to remove the background in such mode. The present observa¬ 
tions of NGC7023 were deemed to have calibration value by 
the Herschel Science Center, and it was decided to obtain an 
additional set of observations at the centre position so that the 
PDR would be observed with a higher S/N. This was achieved 
by including an off-observation in the astronomical observation 
request (AOR), to be performed after three raster positions and 
with zero offset relative to the centre of the map. This resulted in 
the centre of the map being observed a total of 3 times, with the 
redundancy resulting in a higher S/N for that central position, 
which is centred in the PDR. 

The data were processed using version 6.0.3 of the reduc¬ 
tion and analysis package HIPE ( Otij, 2010 1. Each raster posi¬ 
tion, or footprint, consists of 5x5 spectral pixels (called spax- 
els). HIPE produces a data cube for each footprint, which con¬ 
tains a spectrum for each spaxel. Starting from level 1 the cubes 
were further processed, using proprietary tools, to correct for 
minor drifting effects and flux misalignments between scans. At 
this point the cubes are exported into the software PACSman 
dLebouteiller et ali 12012 :1 to measure the line fluxes by fitting 
a Gaussian, and create the final map. The uncertainties in the 
line fluxes are small and on average amount to less than 5% 
for the [CII] and [OI] lines. The relative flux accuracy be¬ 
tween spaxels is lO^fl and for the remaining of the paper we 
have adopted this value or the uncertainty from the fit given 
by PACSman, whichever is higher. Note that we do not de¬ 
tect the [Nil] 122/mr line at any position in the map. The fi¬ 
nal integrated intensity maps of the [O I]63 /mr, [OI] 145 /mi, and 
[CII] 158 /an lines are shown in Fig. 2. 

We have compared ou r meas urements of the [CII] line with 
the fluxes measured by Jobli n et afl (l2010h using the HIF I instru - 
ment, and the ISO/LWS measurements bv lFuente et al J ([2000). 
all compared to the largest beam. The position used for compar¬ 
ison (RA=315.385°, DEC=68.1743°) is located in the clump of 
enhanced emission and is marked at the top right panel of Figure 
2 with a plu s symb ol. All measurements agree within ~10 %. 
lOkada et al.l ( 2013 1 observed the same lines in the NW PDR. 
Their observations consists of a single pointing, while our ob¬ 
servations perform a detail mapping of the region to show the 
morphology of these lines, providing a complete coverage of the 
region and allowing for better calibration. We have compared our 
fluxes to that of lOkada et al] (1201 3l) in the cavity, interface, and 
molecular region. These three positions are marked with squares 
in the top right panel of Figure 2. All values agree within 29%, 
where the [CII] 158 /mi and [OI] 145 /mi fluxes in the interface. 


2 From the PACS spectroscopy performance 
and calibration manual. This can be found at 
http://herschel.esac.esa.int/twiki/bin/view/Public/PacsCalibrationWeb 


and the [O I]63 /mr flux in the molecular region agree within 5%. 
The uncertainty in the lOkada et~ak ( 12013) paper is of ~30%, 
while in our work is typically <10%, and thus both values are 
consistent within the uncertainties. 

In addition to the atomic lines presented in this paper, we 
compare in Section 4 their emission to other relevant cooling 
species of CO, H 2 O, CH + , and C°. In particular, we use the CO 
4<J„<13 transitions in the SPIRE range, and the rotational lines 
of H 2 0-0 S(l), S(2), and S(3) at 17.0, 12.3, and 9.7 /m 1 from 
the Spitzer data. We also use the 360 and 609 /mr lines of C° 
(SPIRE). Since we do not detect any CH + and H 2 0 lines in the 
SPIRE range, we derive an upper-limit to the transitions of these 
speci es that were detected in the Orion Bar (Ber nard-Salas et all 
l2012h . These ob servations were take n from Pilleri et al. ( 201 21) 
for H 2 lines, and lKohler et al] (12014) for the SPIRE lines. 


3. Spatial Distribution 

3. 1. Cooling Lines 

The three top panels in Fig. 2 show the observed distribu¬ 
tions of the [OI]63 and 145/mr, and [CII]158/mr lines. The 
exciting source, the binary system HD 200775 (RA=315.404°, 
DEC=68.1633°), illuminates the PDR from the lower left side 
of the figure (just outside the maps, see also Fig. 3). For each 
panel the beams are shown in the lower left corner and represent 
sizes of 4.5", 8.8", and 11" from left to right respectively. In 
this scale, and adopting a distance of 430 pc to NGC 7023, 10" 
corresponds to a physical scale of 0.02 pc. 

The PDR is detected in the three lines. Their emission arises 
from the surface of the clouds, where the gas is warm (~7 00K- 
70K, iFuente et~ak[ 119991: iFleming et ali 2010t iKohler et~akl 
|20T1 . There is also evidence for condensations of enhanced 
emission within the PDR, where both [OI]63 and 145 /mi line 
maps can resolve two structures (red clumps in the figure). This 
is the first time in NGC 7023 that such structures have been re¬ 
vealed from the emission of the [CII] 158/mi and [OI]63 and 
145 /on lines within the PDR. These structures or knots are di¬ 
rectly illuminated by the star, and coincide with the presence 
of highly-excited CO lines and dust emission (see next sec¬ 
tion). These clumps are better delineated in the [OI] 145 /mi map 
(Figure 2), which reveals two distinct clumps, and a possibly 
third fainter clump. The size of the main structures is around 8" 
(0.015 pc at the distance to NGC 7023). While it is possible that 
these clumps may not be yet be resolved in our observations, we 
find that these clumps cover 9% of the PDR surface area, where 
in this case we consider the PDR area to be any [OI] 145 //m flux 
greater than 0.9x10 7 W tn 2 sr 1 (red, green, and blue colours 
in the upper middle-panel of Fig. 2). 

The PDR is well delineated in the [OI] 145 /mi line 
(Fig. 2). Further away from HD 200775, the emission of 
[O I] 145 /mi drops to zero. The spatial distribution of the 
[O I]63 /mi line is similar, although the emission is more ex¬ 
tended than that of the [OI] 145 /mi line. This is expected as the 
[OI]63/«n line becomes optically thick at lower column den¬ 
sities than the [OI] 145//m line, and once saturated, its emis¬ 
sion will be broader (ICaux et al.Ul999l) . The [CII] 158/mi emis¬ 
sion is even more extended, since it is easily excited. This 
was also seen in a simil a r stud y by our group on the Orion 
Bar (iBernard-Salas et ak[ 120121) . The figure also shows that 
the morphology of the [CII] 158 /mi line is complex and 
structured. Using Herschel/HIFI high-resolution spectra of the 
[CII]158//m line. lJoblin et al.l (1201 Ot) and Berne et al. (2014, in 
prep.) detect different components with a large distribution of ve- 
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Fig. 2. The upper three figures show the observed images for the [O I]63 pm, [OI] 145 /mi. and [CII] 158 /mi lines. The lower panels 
show combination of lines with contour maps over-plotted in white. The bottom images have been convolved to the 158 /mi (largest) 
beam, with the beam size illustrated in the bottom left of each panel. The squares and plus symbols in the upper right panel indicate 
the positions used to compare the [CII] 158 /mi flux to previous measurements in the literature (see Sect. 2). All maps are in flux 
units of 10- 7 W irT 2 sr~'. The eight highest contours are at 100, 91, 82, 74, 65, 56, 47, and 38% of the peak emission. 


locities, indicating that besides the main filaments, other regions 
are contributing to the [CII] emission. 

In the lower three panels of Fig. 2 we plot different com¬ 
binations of lines and their contours convolved to the largest 
beam size (11")- It can be seen that both of the [OI] lines corre¬ 
late spatially very well with each other (lower left panel). The 
middle bottom panel shows that the [CII] emission also de¬ 
lineates the PDR but with a more complex structure than the 
[OI] 145/rm line. The peak emission is slightly displaced and 
broader in the [CII] line than the [OI] 145/mi. This is bet¬ 
ter seen in the last panel where the concentric contours of the 
[OI]63/mi line contrasts with the more irregular emission of 
[CII], including a peak of emission with a heart -shape structure. 

Finally, as mentioned in the previous section, no [N n] line 
at 122 /mi was detected in the PACS map. The 205 /mi line was 
also not detected in the SPIRE maps (iKohler et all 12014 ). Using 
HIFI observations that include the ionisin g sta rs, the Herschel 
Warm and Dense ISM program (WADI. lOssenkopf et all 1201 ll ) 
detects ionised emission only from a shell close to the stars (pri¬ 
vate communication), an area just outside our observed map. 
This could indicate that the gas in the cavity near the PDR is 
mostly neutral and not ionised, and explains why we do not de¬ 
tect the 122/rm line since our map does not include this region 
close to the ionisin g sta rs. Alternatively, if the density is too low 
in the cavity dFuente et all 12000 ) its emission will be too weak 
to detect in our observations. We have derived an upper-limit to 
the flux for this line of l.lxl0~ 7 W m -2 sr _1 . 


3.2. Comparison with CO, H 2 , and PAH emission 


Fig. 3 (left) displays a high-resolution CFHT image of the FU 
1-0 S(l) molecular line at 2.12pm dLemaire et al.Lfl996h . which 
reveals in detail the filamentary structure of the region. Over¬ 
plotted in white and red respectively are contours of the [OI] 
145/rm line and one of the excited rotational 12 CO line (.1=12- 
11) measured in the SPIRE data by iKohler et ali ( 2014 ). We 
note that the images are not convolved as the purpose of this 
panel is to place the emission of [OI] 145/rm line and warm 
CO relative to the filaments. The right panel of the figure shows 
the Spitzer FU 0-0 S(3) 9.7 pm map convolved to the 145/rm 
beam, with white [OI] 145/rm contours. The Hi 0-0 S(3) tran¬ 
sition follows the same distribution as the 1-0 S( 1), but being 
pure rotational line it peaks slightly behind (~1.9") the 1-0 S(l) 
transition. These observations show for the first time that the 
[OI] 145 pm emission in NCG 7023 is displaced with respect to 
the bright filaments that are seen in H 2 . The reason for this is that 
since the H 2 emission closely traces the UV field, its emission 
peaks at the edge of the PDR, while the [OI] 145/r m is m ore 
sensitive to both the gas density and temperature. IKohler et al.l 
(|2014l) have recently mapped the emission of high excitation CO 
lines in the region. It can be seen that this highly excited CO 
emission (left panel in Fig. 3) is located at the same position as 
the [OI] peak. We note that the I2 CO .1=12-11 line peaks at the 
same position as the also observed 17 CO J=7-6 line. Similarly, 
although not shown, these condensations are also detected in 
Herschel/SPIKE and IRAM submillimitre maps that trace dust 
emission (IKohler et all 1201 4 . 


4 






























Bernard-Salas et al.: Cooling lines in NGC 7023 


cn 

CD 

"O 


CJ 

LxJ 

Q 


315.42 315.40 315.38 315.36 315.34 315.32 315.41 315.40 315.3S 315.38 315.37 315.36 

RA (deg) 



68.18 


68.17 


68.16 


68.19 


[Oil 145 contour 


68.180 


68.165 


Fig.3. Left panel: Spatial correlation of the Hj 1-0 S(l) emission map at 2.12 pm ( iLemaire et al.L fl99 6i) compared to the [OI] 
(145/im) and 12 CO .1=12-11 (jKohler et al.L |2014|) contour maps in white and red respectively. To better illustrate the filaments these 
images have not been convolved to the largest beam (in this case, CO). The two orange symbols (squares) mark the position where 
the high excitation CO lines are detected (A), and a nearby position with similar conditions but showing no highly excited CO 
emission (B). These positions are discussed in Sect. 3.2. Right panel: hC 0-0 S(3) map at 9.1pm convolved to the [OI] 145/mi beam 
with its contours in white. In both panels the colour bar intensities are in units of 10~ 7 W m -2 sr 1 , where the four highest contours 
have values of 5.4, 4.7, 3.9, and 3.1 for [OI] 145/on, and of 0.37, 0.3, 0.2, and 0.1 for the CO. 


It is interesting to explore the origin of this highly excited 
CO in relation to the atomic cooling line emission. For this pur¬ 
pose, and to guide the discussion in the next section, we have 
selected two regions: one coinciding with the CO J=12-l 1 peak 
(position A), and a nearby region with no high-J CO emission 
( B ). We note that low-J emission is detected in positions A and B 
( Kohler et all I20l4) . Both positions are labelled in Fig. 3. The 
[OI] 145 pm flux and the [OI] 145 pm/[C II] 158 pm ratio in po¬ 
sition A is respectively 2.2 and 1.7 higher than in position B. 
Although small, these differences are real (larger than the asso¬ 
ciated uncertainties in the [OI] 145/.mi flux and ratio), and in¬ 
dicate variations in the physical conditions between the two re¬ 
gions considered. Such differences can be achieved by changes 
in the radiation field and/or density of a factor ~3 in standard 
PDR grid models (iKaufman et all 120061) . How can such small 
differences explain the presence of high-excited CO in posi¬ 
tion A? The emission of highly excited CO can be explained 
by changes in the gas density if acertain threshold in the condi¬ 
tions is reached. [Kohler et~al] (120141) derived a density of 5 x 10 4 - 
10 6 cm~ 3 and a temperature ranging between 65-130 K at po¬ 
sition A, while at position B they found a lower density (10 4 - 
10 5 cm -3 ) and T=80-120 K. The critical density of the CO lines 
(10 5 —10 7 cm' 3 for .1=4-3 to 13-12 transitions) is higher than the 
[OI] lines (~10 5 cm -3 ). This means that, in the density regime 
between 10 4 6 cm -3 , the CO lines are more sensitive than the 
oxygen lines to an increase in density which could boost the CO 
emission while having little effect on the atomic line emission. 
It is thus possible that differences in densities in the two posi¬ 
tions could thus explain the presence of highly excited CO in the 
clumps. 


Because of the low ionisation potential required to ionised 
carbon (11.3 eV), the [CII] 158 pm line can arise from both the 
neutral (i.e. PDR) and the ionised medium. To validate its use 
as star formation indicator, it is interesting to compare the [CII] 
emission with other known tracers of the star formation process, 
like PAHs dPeeters et all 12004 ). In Fig. 4 we make such compari¬ 
son where we plot the PAH emission traced by the Spitzer/IRAC 
3.6 pm band with the [CII] 158 pm contour map. Despite the 
the fact that the [CII]158/mi line includes emission from dif¬ 
ferent velocity components, it can be seen that the [CII] line 
follows reasonably well the PAH emission, with the peak of 
[CII] 158/mi emission (heart-shaped) coinciding with the PAH 
filament. Whichever regions contribute to the [CII] 158 pm emis¬ 
sion, these must be (mostly) PDR dominated since it traces the 
PAH emission. 

3.3. Line Correlations and Ratios 

Different lines and line ratios have been plotted against each 
other as a function of the distance to the ionising stars HD 20077 
in Fig. 5. Panel a in this figure highlights the different regions 
we consider. These include regions in front of and behind the 
PDR, which in subsequent panels are shown as blue and red 
dots respectively. In the figure the plus symbols are PDR points 
with a flux higher than 75% of the [OI]145//m peak emission 
and are delineated by the black contour. Also shown are se¬ 
lected points from an adopted cut across the PDR starting from 
HD 200775 (labeled numbers). The latter are shown with their 
error bars to illustrate the flux uncertainty in different regions of 
the map. The variation in the flux intensities across the PDR of 
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Fig. 5. Intensity plots for different combinations of line fluxes (in 10 7 W m~ 2 sr _1 ) and line ratios relative to other lines or distance 
to HD 200775. All fluxes have been convolved to the largest beam (11")- The different regions are indicated in the convolved 
[O I] 145 gm map in the first panel and are labelled in the second, where blue represents points in front of the PDR, and red those 
behind it. The peak emission in the PDR is delineated by the black curve in the first panel (those with [OI] 145 pm fluxes higher 
than 75% of the peak intensity of the line), and are plotted as plus symbols in the subsequent panels. The reference points for the 
adopted cut are numbered and shown in the first panel (Sect. 5), and have been plotted as black dots with error bars in subsequent 
panels to give an indication of the uncertainties. 


the [CII] 158 gm, and [OI]63 and 145 gm lines is larger than a 
factor 10, which is higher tha n those observed in the Ori on Bar 
(factor <7) for the same lines (Bernard-Salas et afl l2012l) . 

Panels b and c show respectively the relation between [C n] 
and [OI]63 and 145 pm lines. There is a broad correlation in 
front of and after the PDR in both plots. The likely reason for 
the broad relation is the mixture in the line of sight of the dif¬ 
ferent components that give rise to the [CII] 158gm line (see 
Sect. 3). The relation in panel c is somewhat narrower than in 
panel b. This could be attributed to the fact that the emission of 
the [O I]63 pm line, being optically thicker than the [OI] 145 pm, 
is also extended and follows the [CII] 158 pm emission more 
closely. 

Panel d shows a linear correlation between the [OI]63 and 
145 pm lines behind the PDR, and a more parabolic relationship 
in front of it. In front of the PDR the gas is warmer and op¬ 
tical depth effects become more important. Having an Einstein 
coefficient five times higher, the [OI]63pm line will be more 
self-absorbed than the [OI] 145 pm line, and its overall emission 
will be more extended. Therefore, in that region [OI] 145 pm in¬ 
creases more than the [OI]63 pm line does. 

In the last two panels in the figure (e and f), the 
[OI] 145 pm/[C II] 158 pm and [OI] 145 pm/[01]63 pm ratios are 
plotted respectively against the distance to HD 200775. As we 
approach the PDR the density increases and both line ratios in¬ 
creases rapidly, peaking at around 45". After the PDR, the tem¬ 
perature drops and both line ratios decrease. However, and be¬ 


cause of the multi-component [CII] 158 pm emission, there is a 
broader relation in the ratio. 

4. Cooling budget 

The mid- and far-IR regions are home to the most important 
cooling lines in the PDR. Our Herschel observations, together 
with Spitzer data, have enabled us to estimate the importance of 
these species to the cooling budget. We have thus derived the 
contribution of the [CII] and [OI] cooling lines in the PDR, as 
well as other relevant atomic and molecular lines in the mid- and 
far-IR including; CO, H 2 , and C°. H 2 and CO are the most abun¬ 
dance molecular speciesjn the ISM and contribute also to the 
cooling in PDRs dHollenbach & Tielens . 1999 ). To calculate the 
budget of these species we consider the 0-0 S( 1), S(2), and S(3) 
rotational lines of H 2 at 17.0, 12.3, and 9.7 pm (Spitzer/IRS), 
and CO 4< J„ <13 transitions (SPIRE). In the Orion Bar, C°, 
CH + and H 2 O lines are detected in the PACS and SPIRE spec¬ 
tra and contribute ~ 1 % of the tota l cooling budget (Habart et all 
12010t iBernard-Salas et all l2012h . This is a small factor to the 
total cooling budget in the Orion Bar, but we have explored their 
contribution in NGC 7023 as well. We have detected the 370 
and 609 pm lines of C° in the SPIRE spectrum dKohler et all 
l2014i) . We do not detect CH + and H 2 0 transition in the SPIRE 
spectrum of NGC7023 so have derived upper-limits for their far- 
IR transitions to estimate their contribution. For this comparison 
we choose the position of the [OI] emission peak in the PDR 
(red knots in Fig. 2), and we convolve all the line emission to 
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Fig. 4. Comparison of the PAH morphology as traced by the 
IRAC 3.6pm band intensity in MJy/sr (taken from the Spitzer 
archive), with [CII] 158 pm in white contours. In decreasing or¬ 
der the contour levels are at 10.1, 9.0, 7.8, 6.7, 5.6, 4.5, 3.4, and 
2.2 xl(T 7 W rrT 2 sr -1 . 


the largest beam, in this case the SPIRE CO emission (25"). In 
Table 1 we quote these values together with those of the Orion 
Bar (Bernard-Salas et all 12012 ). which is usually adopted as the 
prototypical PDR. 


Table 1. Cooling line emission as a percentage at the PDR posi¬ 
tion and convolving the emission of all species to a 25 beam. 


NGC 7023 

Orion Bar T 

[C11] 158pm 

19 

8 

[OI] 63 pm 

33 

72 

[OI] 145 pm 

9 

10 

[CII] + [OI] 

61 

90 

Hi 

35 

5 

CO 

4 

5 

C°, HiO, CH+ 

<1 

<1 


' From lBemard-Salas et allll2012h . 


Table 1 shows that in the NW-PDR of NGC 7023 the cool¬ 
ing budget is dominated by the [CII]+[OI] (61%), followed by 
H 2 with 35%, CO (4%), and C° contributing less than 1%. The 
uncertainty in these percentages is around 20% of each value. 
The main contributors to the cooling budget are [O I]63 //m and 
Hi. The overall result is similar to that found in the Orion Bar, 
where the atomic lines dominate the cooling line emission, with 
[OI]63 pm being the strongest. However, the main difference is 
the reduced role of H 2 in the Orion Bar (5%) when compared to 
NGC 7023. As part of the same program, we have studied nine 
other nearby PDRs and find that, like for NGC 7023, H 2 con¬ 
tributes to about 30% (Bernard-Salas et al. in prep.). 

In Figure 6 we show the relative contribution by percent¬ 
age of the [OI]63 and 145 gm, and [CII] 158 pm lines, along 


with that of H 2 . For H 2 we consider the 0-0 rotational transi¬ 
tions S(l), S(2), and S(3) and the vibrational line v=l-0 S(l). All 
maps have been convolved to the [C IT] 158 pm beam. Thus, for 
a given point in the map, the total cooling of these four species 
amounts to 100%. Therefore this figure shows which species 
dominates the cooling in different regions. We ignored the con¬ 
tribution of C°, H 2 O and CH + , but as we can see from Table 
1, their contribution is already negligible in the PDR, where we 
would expect their contribution to be strongest. If we first look at 
each map individually, we see that [O I]63 pm contributes more 
in front of the PDR, and [CII] 158 pm outside it (in front and 
behind). [OI] 145 pm and H 2 contribute more in the PDR, but 
with [OI] 145 pm contribution being higher in the south-east re¬ 
gion (lower-right). The [O I] 145 pm line is more sensitive to the 
gas density than the H 2 lines, and thus contributes more where 
the clumps are located. Reading the percentage from the colour 
bars, we can also compare the relative contribution amongst each 
other. In general, [O I]63 pm is the dominant contributor over the 
entire region (>30%). However, [CII] 158 pm is also an impor¬ 
tant coolant outside of the PDR (contributing above 30%), as is 
H 2 in the PDR (>35%). On the other hand, the contribution of 
[OI] 145 pm does not exceed 12% in the PDR and is lower than 
5% in the rest of the region. This figure delineates the impor¬ 
tance of [O I]63 pm, [CII] 158pm, and H 2 as cooling agents as a 
function of the region. 


5. Modelling 

In this section our goal is to explore whether previously reported 
parameters used to model the dust emission in the PDR (density 
profile, geometrical parameters), and CO lines (gas physical con¬ 
ditions), can also reproduce the [OI] 145 yum emission. This is 
important in our quest to find a model that can self-consistently 
reproduce the dust, molecular, and atomic PDR components. 
Using Herschel photometry and spectroscopy, the dust emission 
and molecul ar content (CO) of the NW PDR has been recently 
modelled bv I Arab et al.l (12012h and iKohler et al j (120141) . In our 
observations the [CII] 158 pm line is a mixture of many compo¬ 
nents, the emission of the [OI]63pm line is optically thick, and 
the [NII] 122pm was not detected. Thus, attempting to model the 
atomic component from our observations will serve to provide 
basic physical conditions of the PDR based on just fitting the 
[OI] 145 pm line. This is why we prefer instead to test whether 
previously and more detailed modelling studies can also explain 
the [OI] 145 pm emission. 

In order to model the emission of the line we have selected a 
cut originating from the ionising stars and intersecting the NW 
PDR at a point of peak dust emission at 250 pm. This is the same 
cut used in dKohler et all 12014 1 for the molecular line and dust 
emission. The observed emission profiles for the three lines, to¬ 
gether with that of the H 2 1-0 S(l) and CO .1=12-11 lines, are 
shown in Fig. 7 (left panel). These profiles are shown at their nat¬ 
ural resolution to maximise the spatial information we can obtain 
from them. It can be seen that [O I]63 and 145 pm, [CII] 158 pm, 
and CO J=12-ll peak at the same position. Having the high¬ 
est spatial resolution, the H 2 profile resolves two filaments (see 
Fig. 3), with the densest filament peaking immediately in front 
of the fine-structure and CO lines relative to the external illumi¬ 
nation. The [CII] 158 pm line has the broadest profile, followed 
by [OI]63pm. 

We have modelled the emission of the [OI] 145 pm line 
with the Meudon PDR radiative transfer c ode (version 1.4.4, 
iLe Petit et alll2006l : lLe Bourlot et afll2012l) . This code assumes 
a plane-parallel slab of gas and dust that is illuminated by an in- 
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Fig. 6. Relative cooling budget of the [OI]63 and 145 fim, [CII] 158 fim lines, and the 0-0 S(l), S(2), S(3) Fb rotational transitions 
and v= 1 -0 S(l) in percentage (%). Adding the contribution of all four species for a given point in the map amounts to 100 %. Other 
species that can contribute to the cooling are ignored in this figure (see Sect. 4). In the first panel, and to guide the eye, a contour of 
the PDR as traced by the [O I]63 fin l emission is shown. The images have been convolved to the largest of the beam sizes of these 
species, in this case the [CII] 158 fim line, shown as a solid circle in the [CII] panel. In addition, and for comparison, we also show 
the largest CO beam size (25") used to calculate the cooling budget in Table 1 as a dashed circle . 


cident radiation field, and resolves iteratively the radiative trans¬ 
fer and thermal and chemical balance to compute the atomic and 
molecular emission in the cloud. As in the papers listed above, 
we have assumed an incident radiation field of G~2,600 on the 
NW PDR at 42” from the star (Choks hi et all 1 19881 


Pilleri et a I 


2012). The RADEX analysis of the CO lines (Kohler et al 


20141) . and the modelling of the dust emission are compatible 
with two maximum gas densities of 10 5 and 10 6 cm~ 3 . Thus, 
for our comparison we consider both possibilities. These studies 
derive a PDR length of 0.2 pc for a gas density 10 5 cm~ 3 and 
0.03 pc for a density of 10 6 cirT 3 , giving both a total column 
density of about 10 22 cm -2 , and we adopt these values in our 
study. We assume the density profile derived from the dust anal¬ 
ysis and adjusted to start at 0.088 pc (42") from HD 200775. The 
density profiles for the two gas densities we consider are shown 
in the middle panel of Fig. 7. We hereafter refer to these mod¬ 
els as model 1 (n=10 5 cm -3 ) and and model 2 (n=10 6 cm -3 ), 
respectively. 


Before we proceed, it is interesting to investigate whether the 
depth of the PDR deduced from the previous studies (lArab et all 
120121: IKohler et all 12014 is compatible with the width and 
length of the PDR we derive from the observations (spatial ex¬ 
tent of the [O I] 145 fim line). In the middle-top panel of Fig. 2, 


the PDR has an approximate width of 20" and a length of 60" 
(as traced by the green colours in the figure), which corresponds 
to 0.04 and 0.12 pc respectively. From model 1 we find that the 
depth is about 2 times the observed length, resulting in a flat ex¬ 
tended geometry. For model 2 the the depth compares very well 
to the projected width resulting in a cylindric geometry. 

The predicted [OI] 145 fim intensities for models 1 and 2 are 
shown in the righthand panel of Fig. 7, where the line emissivi- 
ties have been convolved to the PACS beam at 145 fim. For this 
comparison we consider the conditions of the models in the re¬ 
gion where the predicted unconvolved flux is greater than half 
the maximum flux (48"—52" from the binary system). In this 
region model 1 gives n =10 4 -8xl0 4 cm -3 and T =90-530 K. 
Model 2 gives n =2xl0 4 -2xl0 5 cm -3 and T =90-660 K. In both 
models most of the emission comes from a region with a visual 
extinction Ay = 0.07-1.2. With a 10% uncertainty in the ob¬ 
served line, none of the models reproduces the [OI] 145 fim peak 
intensity. Note that we did not attempt to fit the line, but this 
additional constraint may serve to discriminate between the two 
models. Model 1 predicts an intensity nearly three times higher 
than the observations, while model 2 underestimate the emission 
by just 30%. The maximum gas temperature in both models is 
higher than the gas temperature derived by iKohler et al.l C 20141) 
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Fig. 7. Left panel: Observed non-convolved profile of the atomic and molecular lines as a function of the distance to HD 200775 
along the cut shown in Fig. 5 and described in Sect. 4. Middle panel: Density profiles at two different maximum densities used by 
iKohler et ali ( 2014 1 in their analysis of the molecular emission. Model 1 has a peak gas density of 10 5 and model 2 of 10 6 cm~ 3 . 
Right panel: Observed profile emission of the [O T] 145 pm line compared to the convolved predicted emission by Meudon PDR 
models. 


when modelling the CO J=13-12 to .1=6-5 lines with RADEX 
(140 K). The higher temperature we find is expected since mod¬ 
els predict that the [OI] 145/mi line originates from a slightly 
warmer layer than that of the intermediate J CO excited lines. 

The fact that we find a model that can reproduce the dust 
emission, as well as roughly matching the spatial profile of the 
[OI] 145 pm line is an encouraging step towards a self consis¬ 
tent model that explains both the gas and dust. However, while 
model 2 does not reproduce exactly the peak and width of the 
[OI] emission, we must bear in mind that we are assuming a sim¬ 
ple model with a simple geometry. The H 2 observations reveal a 
more complex structure with several filaments. Also the density 
profile we use reflects a smooth PDR, while our own observa¬ 
tions clearly reveal the presence of dense structure inside of the 
PDR. Our cut is also arbitrary. These complex structures will af¬ 
fect the [OI] emission profile (width, peak position). Given these 
caveats, it is therefore unsurprising that we do not have a perfect 
fit. In fact, in order to explain the intensity peak of the excited 
CO emission, a steeper density gradient at the PDR edge over a 
small scale is needed, to have both a large amount of warm and 
dense gas. The CO lines can be reproduced using a high pres¬ 
sure model (P =10 8 K cm" 3 ) and relatively small PDR width 
and depth of 2-3" (Joblin et al. in prep.). Such a model can also 
predict the maximum peak intensity of the [OI] 145/mi line. A 
lower pressure is however needed to reproduce the full extent of 
the PDR that we have presented in Fig.2. 

Finally, and since the Meudon code does not takes into ac¬ 
count the ionised region, we have modelled the emission of the 
[Nil] 122 and 205 pm lines using the radiative transfer code 
Cloudy dFerland et all 119981 1. We use the same density profile, 
peak density, and PDR length as in motleQ 2. Fo r the radia- 
tion field, we assume a Kurucz spectrum ( Kuruc z] 1 19931) with 
Teff =15,000 K corresponding to the spectral type ( B3V-B5V) 
of the illuminating binary star ( lAlecian et all 12008 .) but atten¬ 
uated in order to have an incident radiation field of Go~2,600 


3 Using the parameters of model 1 does not change the conclusion. 


at the edge of the PDR dPilleri et al.. 2012 ). While there is a 
cavity, the ionising sources are young and likely to have some 
material around the star. This attenuation of the radiation field 
is likely caused by the envelope of t he star and the dilut i on ef¬ 
fect. We also adopt ISM abundances (Savage & Sembach, 1996; 
iMever et al.L1 1 99711 1 9981) . Between the star and the PDR there is 
a cavity. The density profile is adjusted to start at 0.023 pc (11") 
from HD 200775. The intensity of the nitrogen lines at 122 and 
205 pm are sensitive to the electron density in the ionised region, 
and on the incident radiation field, but we do not detect either 
line. The upper limits we have derived for these lines (1.1 x 10 7 
and 10 9 W nV 2 sr~* respectively) are compatible with two sce¬ 
narios for the cavity, one with an ionised region and a density 
H<100 cm -3 , and another one with no ionised region and a den¬ 
sity of 500 cm -3 . 

6. Summary and Conclusions 

We have presented Herschel /PACS spatially resolved observa¬ 
tions of the [CII] 158 pm, [OI] 63pm and 145 pm lines of the 
NW PDR in NGC 7023. This has enabled us to study the emis¬ 
sion of these cooling lines in relation to the morphology of the 
region. We summarise here the major findings; 

1. The emission of the atomic cooling lines trace the cloud 
surface that is directly illuminated by the binary system 
HD 200775. Here the gas is warmer, and these lines are as¬ 
sociated with the filamentary structure revealed in higher- 
resolution H 2 and PAH maps. In the PDR, the peak of cool¬ 
ing line emission corresponds to the presence of condensa¬ 
tions that are seen from high level rotational CO lines, and 
Herschel and submillimetre emission of dust. 

2. By comparing the role of different coolants we find that the 
[OI]63 and 145 pm and [CII] 158 pm lines account for 61% 
of the cooling in the PDR, with [OI]63 /mi contributing to 
33% of the emission. H 2 also contributes significantly with 
35%, CO with 4%, and other atoms and molecules (C°, H 2 O, 
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CH + ) contributing less than 1 percent. Looking at the rela¬ 
tive contribution of the main cooling agents over the region, 
we find that while the [OI]63 pm dominates the cooling, H 2 
also contributes significantly in the PDR, and [CII] 158 pm is 
also important in front of and behind the PDR. This high¬ 
lights the importance of [OI]63pm and PL as a coolant in 
the PDR, and of [CII] 158pm in lower excited regions. It is 
interesting to notice that while in the Orion Bar the atomic 
lines dominate the cooling, PL contributes to only 5%. 

3. The [O I]63 and 145 pm maps spatially resolve these conden¬ 
sation into two distinct structures. Furthermore, HerscheF s 
high-angular resolution shows that the [OI] 145 pm emission 
peaks slightly farther away from the ionising stars than the 
PL emission which traces the edge of the PDR. This is ex¬ 
pected since [OI] 145 pm is more sensitive to the gas density. 

4. [OI]63 and 145pm and [CII] 158pm peak at the same po¬ 
sition. The [CII] 158pm line is more extended than the 
[OI]63 and 145 pm lines because [CII] 158pm is easily ex¬ 
cited and presents a more complex structure. Thi s is c onsis- 
tent with HIFI observations of the region by dJoblin et all 
120101) which reveal different velocity components for the 
[CII] 158pm emission, indicating that different regions are 
contributing to the emission of this line. 

5. The emission of highly excited 12 CO .1=12-11 correlates with 
the peak of emission of the cooling lines. We find that emis¬ 
sion from the atomic cooling lines and line ratios at this po¬ 
sition differ by a factor of two compared with the conditions 
of a nearby position with no CO J=12-11 emission. This dif¬ 
ference in the cooling line emission indicates a change in 
the physical conditions (density varying from 10 4 to 1 0 5 6 
cm -3 ). In this density range, the CO lines are more sensitive 
than the atomic lines and could explain the presence of warm 
CO in this region. 

6. Using a density profile derived to reproduce the dust emis¬ 
sion and physical parameters from the analysis of CO lines, 
we have modelled the emission of the [OI] 145 pm line us¬ 
ing two different peak gas densities using the Meudon code. 
The best model predicts the [OI] 145 pm emission to within 
30%, of the intensity at the peak of the emission, occurring at 
Ay =0.07-1.2. In this region the conditions are n = 1.7xl0 4 - 
1.8xl0 5 cm- 3 , T =90-660 K. The PDR depth along the line 
of sight needed to reproduce the absolute [OI] 145 pm inten¬ 
sity is 0.03 pc, which is comparable to the observed width of 
the PDR in our map. 

7. We did not detect the ionised line of [Nil] at 122pm. Our 
Cloudy model indicate that the ionised shell may have a low 
density, in which case, the [N 11 ] 122 pm line is too weak to be 
detected in our observations. Upper-limits to this line and the 
[N 11 ] at 205 pm line are consistent with two scenarios for the 
cavity, one with an ionised region and a density n < 100 cm -3 , 
and one with no ionised region and a density of 500 cm -3 . 
This is supported by supporting HIFI observations which re¬ 
veal that the [N 11 ] at 205 pm line is only emitting in a shell 
close to the star, and indicates the cavity between the star and 
the PDR is mostly filled with non-ionised gas. 

The morphology of the [OI]63 and 145 pm and 
[CII] 158 pm lines give an unprecedented view into the 
variation of the physical conditions, energetics, and cooling 
across the NW PDR. The clumps of emission we observe 
provide a means to understand how stellar radiation interacts 
with the gas, and in the future it will be interesting to model 
whether such clumps will be photo-evaporated or form stars. 
Given the complexity of the region, and that we have assumed 


a simple planar model to reproduce the [OI] 145 pm using 
parameters derived from modelling the dust and CO emission, 
this study offers an encouraging step forward in producing a 
self-consistent model that can explain the atomic, molecular, 
and dust emission. Moreover, spatially resolved studies like 
this can help us calibrate the use of the [C II] 158pm and [OI] 
lines to trace star formation dStacev et all 1991 1 pMei ierink et all 
120071), and understan d the origin of the so called [CII] deficit 
( Luhman et all 120031) in most luminous galaxies. Because of 
its brightness, the [CII] 158pm line provides an unmatched 
opportunity for redshift determinations and source diagnostics 
of far distant systems, especially with ALMA. 
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